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We propose a scheme to simulate one-dimensional XXZ-type Heisenberg spin models with compet- 
ing interactions between nearest-neighbors (NNs) and next-NNs in photon-coupled micro-cavities. 
Our scheme, for the first time, exploits the rich resources and flexible controls available in such a sys- 
tem to realize arbitrarily adjustable ratios between the effective NN and next-NN coupling strengths. 
Such powerful capability allows us to simulate frustration phenomena and disorder behaviors in f-d 
£S| , systems arising from next-NN interactions, a large class of problems of great importance in con- 

densed matter physics. Our scheme is robust due to the lack of atomic excitations which suppresses 
spontaneous emission and cavity decay strongly. 
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Quantum simulation is an important application of quantum information science. As a promising physical system 
t-H ■ for quantum simulation, ultracold atoms trapped in an optical lattice [l|, 0] have many appealing properties such as 
long coherence times, the possibility of simultaneous initialization of a large number of particles, and tunable effective 
coupling strength within a large range. However, it is difficult to focus a laser beam on a single atom due to the 
short optical lattice period which makes it challenging to realize single-site operations [3| . To address this challenge, 
recently arrays of photon-coupled microcavities have been suggested as an alternative quantum simulator P4l9| . An 
artificial system on a micro-chip, microcavities can be fabricated in desired structures and dimensions with great 
precision, a significant advantage that allows to realize both single-site operations and neighboring-sites interactions 
easily Q. 

In previous research on quantum simulation, the focus has been on realizing various spin models with on-site and NN 
i ^ i . interactions, the most widely used model in studies of condensed-matter physics problems such as quantum magnetism 
[20L I2II ] . For instance, it was demonstrated that an array of coupled cavities with one atom in each cavity can be used 
t-H \ to simulate the anisotropic (XYZ) Heisenberg spin-i model. It was further shown that the XXZ-type Heisenberg spin 
model of any high spin can be realized using a cavity array with a number of atoms in each cavity fal. The XXZ-type 
Heisenberg spin model can also be realized when the spin states are represented by polaritons flOlllTI . Much of the 
recent research was reviewed in Q . However, it is well known in condensed- matter physics that many important physics 
jy-j ' and exotic phenomena arise due to the long-range nature of interactions between spins. A particularly important 
model to capture this is a spin model with both (and possibly competing) NN and next-NN interactions, which is 
known to reveal many important physics such as Tomonaga-Luttinger liquid states and spin-Peierls states |2lH23| . It 
is even suggested that it could be used to study the mechanism for iron-based superconductors [24j . The importance of 
I ■ such frustrated spin systems has also been recognized in the context of quantum information. For instance, behavior 
of quantum entanglement in frustrated spin systems was studied by a number of authors (25l - [29| and interesting 
and deep connections between entanglement and phase transitions were dicovered in these systems. Unfortunately, 
physical implementation for a spin-chain quantum simulator with both tunable NN and next-NN interactions has not 
been available, mainly due to the technical difficulty in realizing controllable next-NN interactions. For instance, in 
ultracold atomic systems the adjustable spin interactions are usually realized by controlling the wave function overlap 
between neighboring sites. This technique is not useful in engineering next-NN interactions since the wavefunction 
overlap falls off exponentially with distance. 

In this work, we show how one can implement an effective spin model in a microcavity system with both NN 
and next-NN interactions. By taking advantage of the many controls available and using a smart idea of interaction 
cancellation and enhancement, we can adjust at will the ratio between the NN and next-NN interaction strengths. This 
then allows for the first time to simulate a large class of condensed-matter physics problems in which the competition 
between NN and next-NN interactions plays an essential role. Moreover, in our approach, cavity field is always kept 
in ground state and atoms remain in two long-lived states. Since excited atomic levels and cavity photon modes only 
appear in virtual processes, spontaneous emission of internal states and cavity decay are strongly suppressed. 

We consider an array of cavities that are coupled via exchange of photons with one 4-level atom in each cavity 
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FIG. 1: (color online), (a) A one-dimensional array of coupled cavities with one 4-level atom in each cavity, (b) Involved 
atomic levels and transitions. 



(Fig.l). Such a model can be realized in several kinds of physical systems such as microtoroidal cavity arrays [l7^ . 
photonic crystal defects 18] and superconducting stripline resonators [19j. Two long-lived levels, |1) and |2), represent 
the two spin states | I) and | f) for the effective spin. Together with two excited states |3) and |4), they form two 
independent A level structures. We denote the subsystem consisting of levels |1), |2), and |3) A , and that consisting 
of |1), |2), and |4) A*. 

In previous research, a simpler atomic level configuration with one A structure was used to realize effect spin models 
with NN XXZ interactions [8, 9]. In these models, interaction strength between spins decreases rapidly with distances 
and next-NN interactions are negligibly small compared to NN interactions. In our study, by cleverly adjusting 
relevant experimental conditions, we can make interactions arising due to the two independent level structures A a 
and Af, add up or cancel each other depending on the phases of the control lasers. This is the key idea that allows us 
to realize arbitrary ratios between NN and next-NN interaction strengths. 

We now derive the effective Hamiltonian of the photon coupled microcavity system in Fig. 1 (a) when each atom 
couples with two cavity modes and four external lasers (Fig. 1(b)). Suppose two cavity modes, uj a and u>b, are close 
in energy to the transitions |1) <H- |3) and |2) ■<-» |4), and thus drive these two transitions with strengths g a and gi, 
respectively. Further, we apply four lasers with frequencies u>\, u>2, W3 and LU4 to drive the transitions |1) o |3), 
|2) <-» |4), |2) O 1 3) and |1) -h- |4), with Rabi frequencies Sli, ffee - " 1 ^ , 3) Q 4: e~ lnj (the index j represents the jth 
cavity), respectively (without loss of generality, here, assume all the fii are real). Here, we have modulated the phases 
of lasers u>2 and W4 on purpose, and as will be seen later this phase modulation plays a key role in our scheme. All 
the transitions are assumed to be large detuned, i.e., the magnitudes of the detunings £31 = OJ31 — u> a , 642 = ^42 — Wb, 
A31 = W31 — uji, A42 = W42 — u>2, A32 = W32 — 0J3, and A41 = W41 — W4 ( u>^ v is the energy difference between level 
\jjb) and \v)) are much greater than the transition strengths |g a |, \gb\ and = 1,2,3,4). Under these conditions, 

in the rotating frame, we can perform a standard adiabatic elimination of the atomic excited states |3) and |4) (30j 
and obtain an effective Hamiltonian 
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where 8\ = — A31 = 5^ — A42, 5 2 = £31 — A32 = £42 — A41, and S3 = Si — S 2 are the two-photon detunings, 

M = + £), a 2 = a^(£ + £), 4. = + £), Si = + 52 = ^(ait + 

and -B3 = 2 (zr" are * ne effective coupling strengths, J a , J b are the tunnelling rate of photons between 

neighboring cavities(all assumed to be real). Notice that the effective interaction strengths A\ and A2 are induced 
by the level structure A , whereas B\ and B2 arise due to A b , and there is a phase factor of (— 1) J from the Rabi 
frequencies f^e -17 ^ and f^e -1 *^' in the jth cavity. 

We can now derive the effective spin interactions using the virtual-transition induced by effective Hamiltonian in 
Eq. (TTJ). To avoid excitations of real photons and ensure that all two-photon transitions are independent when we 
derive the effective spin Hamiltonian we further require that all two-photon transitions be large detuned, i.e., 

\Si\, |*3 - S 2 \ » \Ai\, \B t \, \J a \, \J b \, ||J|(i=l,2,3). We define the spin operators 5| = i(|2 J )(2,| - 

Sj~ = |2j-)(lj-|, SJ — \lj)(2j\. Before proceeding, we simplify the Hamiltonian in Eq. (fTJ). First, when we consider the 

third term up to second order, it is J^j A\)S^ so the third term and the fourth term only contribute to effective 

local magnetic field. We temporarily neglect them. Second, the fifth term is a perturbation term that modifies the 

2222 
detunings of the two- photon transitions, i.e., 5 a \ — Si + j 2 -, S a2 = #2 + j^i *6i = *i + j^i *&2 = S 2 + Jj-- Finally, we 
assume periodic boundary conditions and take advantage of Fourier transformation to diagonalize the photon coupling 
terms. Defining djibj) = ^= X^feLi Fj k c k (d k ), where N is the total number of the cavities, Fj k = exp(-i^jk) and 

J2f=i FjkFji — N5~u (Ski is the Kronecker function), we have 

E J a{a]a 3 +i + Oj-at +1 ) = ^T afe cj,c fe , (2) 



Y, Jb{b]b J+1 + b s b] +1 ) = T bk d\d k , (3) 

j k 

where T ak = J a £\j j^{F* k F j+ ij+F ok F* +ll ) = 2J a cos{^-k) and T bk = 2J b cos(^k). The Hamiltonian in the rotating 
frame reads 

H = - J2i[Mhi - S*y 5alt + A 2 Sfe iS ^}%c k e- iT °* + H.c.} 

(4) 

" ^{(-tyiBih 1 ! + S])e iSnt + B 2 Sje^ t }^Ld k e- lTb » + H.c.}. 



When the large detuning conditions \S ai \, \S bi \, \S a \ — 8 a2 \, \8 b \ — S b2 \ ^> \Bj\, \ J a \, \ J b \ are satisfied, all terms are 
independent when we adiabatically eliminate the photon states, and we obtain the effective spin Hamiltonian 
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It is clear from Eq. (5) that effective spin interactions arise between any pair of lattice sites j and I due to nonlocal 
modes Ck and dk- However, due to the large detuning conditions |<5 M i| 3> [T^k] = |2J^cos(^fc)| (/i = a, &; i = 1,2), 
the interaction strength drops quickly with the site distance \i — j\. To see this, we expand the following term 
for interaction strengths in Eq. (5) to second order in 7j 1 fc/<W Making use of the relations J^k-FjkF^ — NSji, 
J2k F jk Ft k ■ 2co S (f k) = N(6j, l+1 + 8 J+1 ,i), Ek FokFtk • (2cos(f k)f = N(2~6 M + 5 u+2 + S j+2 ,i), we obtain 
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Substituting the expanded terms in Eq. (5), we obtain the following effective spin Hamiltonian to second order in 

^ °al S bl 



J 



[#^ + (-l) 1 #5 2 2 ](5+5r +1 + 5r5t +1 ) 
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J 2 J 2 (7) 
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+ [#^2 + (-l) 2 #S 2 2 ](^^+2 + ^/ +2 ) 
°a2 °62 

+ hjS*}, 

where hj is the strength of local effective magnetic held 31]. It can be see from Eq. (7) that the effective spin 
interaction strength terms proportional to At and Bi arise from the atom's level structures A a and Af, separately. 
For both A and A&, the next-NN interaction strength is 1 order of magnitude (in J/S) smaller than that of NN 
interaction. Therefore, if only one of the level structures A a and Af, was present, it can be seen from the above 
Hamiltonian that the next-NN interaction would have been negligibly weak compared to NN interaction because their 
ratio is J^/S^i <C 1. For this reason, the next-NN interactions are omitted in previous work [H, H|. However, in 
our system both level structures A a and Af, contribute to the effective spin interactions. When calculating the total 
effective interaction strength, we see that contributions due to A a and Af, tend to cancel for NN interactions but 
add up for next-NN interactions, due to our careful choice of the control laser phases. By using this smart idea of 
interaction cancellation and enhancement, we can arbitrarily adjust the ratio between NN and next-NN interaction 
strengths. Finally, we simplify the above effective spin Hamiltonian and write it in the following form 



H 



E 



(8) 



Here, the effective coupling strengths are J x = 2(jtA 2 2 - 4f-Bf}, J 2 = 1{4tAl + 4j-B%), Xi = 2{j^A\ - jfB 2 ), 

al hi o.1 b2 al fal 

A2 = 2(fA 2 + fB 2 ). 

The interaction cancellation and enhancement is clearly reflected in the above expressions for J\ and J 2 , effective 
interaction strengths between NN and next-NN spins in the Hamiltonian Eq. ©. Obviously, by choosing close values 
-A 2 and 4r-B 2 (i = 1, 2), we can make J 2 comparable to or even greater than J±. We can also adjust the relative 
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signs of J\ and J 2 easily. 
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FIG. 2: Atoms in cavities interact via exchange of virtual photons of cavity modes uu a ((a)) and uut ((b)). The two spin states 
| 4-) and | t) are represented by two long-lived atomic levels |1) and |2). Ai, Bt(i — 1,2) are the effective coupling strengths, 
5ai, Su are the detunings, J a , Jb are the tunnelling rates of photons between neighboring cavities. 



To gain a better intuitive understanding of the critical role that interaction cancellation and enhancement plays 
in the physics underlying the effective Hamiltonian in Eq. ©, we present a simple physical picture in Fig. 2. Fig. 
2 (a) describes interactions induced by the two-photon transitions of the A a level structure. In this process, atoms 
interact via the exchange of virtual photons of the cavity mode uj a . Similarly, Fig. 2(b) describes interactions induced 
by the two-photon transitions of the A;, level structure, where atoms interact via the exchange of virtual photons of 
the cavity mode 0Jb- The difference is, due to the phase modulation of the control lasers, in Fig. 2(b), the amplitude 
of emitting or absorbing a photon is site-dependent, i.e., (— l) J+1 i? i (i = 1, 2) in the jth cavity. 

Now, we consider the first case in Fig. 2(a), which generates ZZ interactions. The probability amplitude that the 
first atom emits a virtual photon with frequency u a is -j^, where S a i is the detuning. Then the photon tunnels 
to the second cavity with amplitude j*-. If it is absorbed with amplitude (— Ai), the strength of NN interaction is 
-yY"/ a j-(~ A-i) — jr^Al > 0. The photon can further tunnel to the third cavity with probability amplitude j^, and be 

absorbed with amplitude (—Ax). This results in a next-NN interaction with strength ^-l) = > 0. 

In Fig. 2(a), similar processes occur and they can be understood by replacing A\ and b a \ with A^ and b a i in the above 
discussion. It generates XX interactions. In Fig. 2(b), due to the site-dependent modulation of the phases of the laser 
u)2 and u>4, the NN sites have different amplitudes, Bi and {—Bi) (i = 1,2). So the effective couplings of the NN 

and next-NN interactions are ^^-(-5;) = -jr-Bf < and j^j^j^B, = ^-Bf > 0, respectively. All the effective 

couplings have a factor of 2 because the atoms that emit and absorb the virtual photon can be switched. We can 
clearly see that the total interactions of the next-NN sites enhance while the total interactions of the NN sites cancel, 
leading to the effective Hamiltonian in Eq. ([5]). 

It should be noticed that the effective Hamiltonian in Eq. ([8]) is highly tunable, because it is determined by 6 free 
parameters (2 parameters on frequency detuning and 4 parameters on intensity of control laser). As long as the large 
detuning conditions are satisfied, the interaction strengths in effective Hamiltonian Eq. ([5]) can be adjusted at will 
allowing to realize arbitrary ratios between J\j Ji and A1/A2. With its large realizable parameter space, it is then 
convenient to use this model to simulate a large class of XXZ Heisenberg spin chain problems in which competing 
interactions between NN and next-NN play an essential role. 

As an example, we discuss how our system can be used to simulate frustration phenomena in Id condensed matter 
physics. By introducing a new stark shift @], the effective local magnetic field can be set to 0. When we choose system 
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FIG. 3: (color online). Time evolution of occupation probability p(lj) of state \lj) (j = 1,2,3,4), calculated using the 
full Hamiltonian (solid blue line) and effective spin Hamiltonian (dashed red line), for parameters Ai = A? = 0.1 GHz, 
Bi = B 2 = 0.096 GHz, A 3 = B s = 0.02 GHz, 5i = 4 GHz, S 2 = 3 GHz, <5 3 = 1 GHz, gl/5 3 i = gl/S i2 = 0.1 GHz, and 
J a = Jb = 0.2 GHz. 

parameters such that J2 = A2 > and J\ > 0, Aj > 0, the Hamiltonian if is the well-known anisotropic Heisenberg 
model with competing interactions ( "frustration" ) originally studied by Haldane [22( . Haldane predicted four phases 
for this system at zero temperature: spin fluid phase, Neel phase, dimer phase, and bound phase. When the ratio of 
J2/J1 is increased, the interesting phenomenon of spontaneous dimerization takes place. 

In coupled cavities systems, addressing of individual cavities is available. Therefore, one can also simulate a 
frustrated spin chain with disorder. Yusuf et. al. predicted the random singlet phase and the large spin phase in 
random antiferromagnetic spin-1/2 chains with NN and next-NN couplings. A strong next-NN coupling will drive the 
system to a large spin phase [32l |. 

Aside from simulating frustrated spin problems of interests in the conventional context of condensed matter physics, 
our system with both tunable NN and next-NN interactions is also very valuable for quantum information studies. 
Spin chains with NN interactions are known to be useful for various quantum information tasks. For instance, 
they can be used to transfer quantum states [33| . Whether a spin chain with both NN and next-NN interactions 
can transfer quantum states more efficiently is an interesting topic to investigate, and our system provides a physical 
implementation for such a system readily. Perhaps more intriguingly, our system offers a powerful quantum simulation 
tool to study the deep connections between quantum entanglement and phase transitions in frustrated spin systems. 
This important subject has been discussed theoretically [26-29], however no physical implementation of such a system 
with controllable NN and next-NN interaction strengths was available to check the findings in these preliminary 
theoretical research. Our system for the first time makes it possible to physically simulate relevant frustrated spin 
models. Because of the large realizable parameter space and ability to address individual spins in our system, we 
can simulate a large class of phase transitions and also study detailed information on entanglement in the system. 
This then allows to verify previous results on entanglement and phase transitions in frustrated spin systems and also 
simulate more complicated problems beyond the reach of conventional numerical methods. 

Our derivation of the effective spin Hamiltonian in Eq. (J5J involves a large number of assumptions and approxi- 
mations. To check the validity of these approximations and the accuracy of the results, we numerically simulate the 
dynamics generated by the full Hamiltonian in Eq. (1) and the effective model described by Eq. © and compare the 
results. 

As an example, we consider four atoms in four cavities, initially in the state |li) ® |2a) ® 1 23) <g) 1 24) corresponding 
to a spin state where only one spin points down, | 4-ttt) ■ We calculate the time evolution of the probability of atom 
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j in the |1) state, i.e., the occupation probability p(lj) of state (j = 1,2,3,4) using the full Hamiltonian in 
Eq. (1). This corresponds to the probability of spin j pointing down, p(ij), which we calculate using the effective 
spin Hamiltonian in Eq. (jSJ. Fig. 3 shows p{lj) and p(jj) for an effective Hamiltonian® with Ji = 0. 0326MHz, 
J 2 = 0.0516MHz, Ai = 0.0187MHz and A 2 = 0.0223MHz |34[. To demonstrate the high tunability of our model, we 
have deliberately chosen the system parameters such that next-NN interactions are stronger than NN interactions. 
Due to periodic boundary conditions, the results of p(l2) and p(li) are the same. As shown in Fig. 3, the numerical 
results of the Hamiltonian Eq. (1) and those of the effective model agree with each other reasonably well. Therefore, 
our effective model is valid. 

Quantum simulation based on our effective spin model is adversely affected by nonidealities in the system, most 
notably the cavity decay and limited lifetime of the excited atomic states. To reliably observe exchange dynamics 
due to next-NN interactions, the coupling strengths J a and X a should be at least one or two orders of magnitude 
larger than the cavity decay rates and the lifetimes of the atomic upper levels. Defining Q = max(51i, 0,2, ^3, ^4) 
, g = max(g , <#,), J = max(J , J&), A = min(A3i, A32, A41, A42), and S = min(5i, #2), one can see that the mean 
population of the atomic excite state |3) (or |4)) can be approximately written as |fi/A| 2 . This results in an effective 
decay rate |f2/A| 2 r£ with Te the linewidth of the upper level. Similarly, the effective cavity field decay rate can 
be expressed as |f2<7/(A<5)| 2 rc, where Tc describes the cavity decay of photons [35J. Since the coupling J a and X a 
is approximately 4J 2 fi 2 # 2 /(A 2 6 3 ), this leads to the constraints T E « ig 2 J 2 /S 3 and T c « -iJ 2 /S. Without loss 
of generality, we assume 5 > 2 J with 2 J the bandwidth of the photons in the cavity array. Thus, to ensure that 
the photons tunnel between cavities before they decay, the condition Tc « J must be satisfied. Additionally, for 
<< 4g 2 J 2 /5 3 , high cooperativity factors g 2 / (TcT E ) and high ratios of coupling to dissipation g/T E are favorable. 

Experimentally, the requirements on parameters discussed above can be fulfilled in microcavities with high quality 
factors. For toroidal microcavities in [17j . the predicted critical atom number can approach 10 -7 , which results in 
a cooperativity factor g 2 /(TcT E ) ~ 10 7 and a ratio g/T E ~ 10 3 . For a strongly coupled single quantum dot-cavity 
system g ~ 20GHz can be achieved with T c = 24GHz and T E = 18MHz, for which we have g 2 /(T c T E ) ~ 10 3 
and g/T E ~ 10 3 . In both systems, the photons can transfer between different cavities via their evanescent fields 
or optical fibres. Thus, our effective spin model with next-NN interactions can be readily and reliably realized in 
presently available systems. 

In this work, we have shown that one dimensional frustrated spin models can be simulated using atoms in photon 
coupled cavities. By choosing the phases of control lasers appropriately, it is possible to make interactions arising 
from different level structures add up or cancel each other, which then allows to adjust the ratio between NN and 
next-NN interaction strengths at will. Quantum simulation using our system is robust and reliable since the atoms 
are only virtually excited and no real photon absorption and generation are involved and thus dissipative processes 
such as photon decay and spontaneous emission are strongly suppressed. 

This work was funded by National Fundamental Research Program (2006CB921900), the Innovation funds from 
Chinese Academy of Sciences, and National Natural Science Foundation of China (Grant No. 60621064, No. 10874170, 
No. 10875110, and No. 60836001). Z.-W. Zhou gratefully acknowledges the support of K.C.Wong Education Foun- 
dation, Hong Kong. 
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